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Abstract. Let D C -R'' be a bounded domain and let 

L = iv-aV + & - V 

be a second order elliptic operator on D. Let be a probability measure 
on D. Denote by £ the differential operator whose domain is specified 
by the following non-local boundary condition: 

Vc = {feC''(D): f fdi. = f\aD}, 

J D 

and which coincides with L on its domain. Clearly is an eigenvalue 
for £, with the corresponding eigenfunction being constant. It is known 
that £ possesses an infinite sequence of eigenvalues, and that with the 
exception of the zero eigenvalue, all eigenvalues have negative real part. 
Define the spectral gap of jC, indexed by v, by 

7i(i^) = sup{Re A : 7^ A is an eigenvalue for £}. 

In this paper we investigate the eigenvalues of C in general and the 
spectral gap in particular. 

The operator £ and its spectral gap 71 (i^) have probabilistic signif- 
icance. The operator £ is the generator of a diffusion process with 
random jumps from the boundary, and 71(1^) measures the exponential 
rate of convergence of this process to its invariant measure. 



1. Introduction and Statement of Results 
Let D C R'^ be a bounded domain and let 



L = -V -aV + b-\/ 
2 

1 
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be a second order elliptic operator on D. We will assume that a = {aij}fj^i 
is positive definite with entries in C^'"(M'^) and that b = (61,..., 6^) has 
entries in C^'"(M'^), for some a € (0,1]. Note that we have written the 
principal part of the operator L in divergence form. This has been done for 
convenience and, in light of the above conditions on the coefficients, without 
loss of generality. We will assume either that D has a C^'°-boundary or that 
D = Di X ■ ■ ■ X Dk, and L = "^^^i Li, where Lj is defined on Di and Di 
has a C^'"-boundary. This latter situation allows in particular for the case 
of on a cube. 

Let 1/ be a probability measure on D. Denote by C the differential operator 
whose domain is specified by a non-local boundary condition as follows: 

Vc = {f€ C\D) : f fdv = f\9r>}. 
Jd 

and which coincides with L on its domain. (Non-local boundary conditions 
in the spirit of the one above in the context of parabolic operators can be 
found in the physics literature on "well-stirred" liquids. See [11] and [5].) 

Clearly is an eigenvalue for C, with the corresponding eigenfunction 
being constant. It is known that C possesses an infinite sequence of eigen- 
values, and that with the exception of the zero eigenvalue, all eigenvalues 
have negative real part (see Theorem BP below). Note that the operator 
C depends on the measure u through its domain of definition. Define the 
spectral gap of C, indexed by z^, by 

(1.1) 71 ('^) = sup{Re A : / A is an eigenvalue for C}. 

In this paper we investigate the eigenvalues of £. in general and the spectral 
gap 71 (z^) in particular. The operator £ and its spectral gap 71 (z^) have 
probabilistic significance which we now point out. 

Let G^{x,y) denote the Green's function for L, defined by 

/•oo 

G^{x,y)= / p^(t,x,y)dt, 

JO 
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where p^{t, x, y) is the Dirichlet heat kernel for L — in D, oi equivalently, 
as a function of y, p^{t, x, y) is the transition subprobability density for the 
diffusion process Y{t) in D corresponding to L, starting from x ^ D and 
killed upon exiting D. It was shown in [Ij that there exists a Markov pro- 
cess X{t) in D which coincides with the diffusion Y{t) governed by L until 
it exits D, at which time it jumps to a point in the domain according to 
the distribution u and starts the diffusion afresh. This same mechanism is 
repeated independently each time the process reaches the boundary. This 
process is called a diffusion with random jumps from the boundary. In light of 
the above probabilistic connection, from now on we will refer to the measure 
V appearing in the definition of C as the jump measure. Denote expected 
values corresponding to this process starting from x ^ D hy E^. Let V{D) 
denote the space of probability measures on D. Under the smoothness con- 
ditions stated above, the following theorem was proven in Theorem 1 and 
the remark following it]. 

Theorem (BP). Let X he the diffusion with random jumps from the bound- 
ary corresponding to L and v . 

i. There exists a unique invariant measure /i for the process. It has a den- 
sity, also denoted by /i, which is given by 

^^^^ ^ J^G^{x,y)du{x) 

The map Inv: 'P{D) V{D), defined by Inv{u) = /i, is continuous in the 
topology of weak convergence of probability measures. 

ii. The operator C possesses an infinite sequence of eigenvalues. Further- 
more, 

lim ylog sup ||E,/(X(t)) - / /d^lU = < 0, 

t^oo t feL--{D), ||/||oc<l Jd 

where 71(1^), defined in (jl.ip . is the spectral gap of C 
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Remark. Actually, part (ii) of Theorem BP was proved in p] for a more 
general problem, where the jump measure from the boundary is allowed to 
depend on the boundary location. 

We now turn to the analysis of the eigenvalues of L in general and of 
the spectral gap of L in particular. Note that by Theorem BP, the larger 
the spectral gap, the faster is the rate of convergence to equilibrium for the 
diffusion with random jumps. 

We begin with a very special case of jump measure v where the eigen- 
values (and eigenfunctions) of L can be completely characterized in terms 
of those of L with the Dirichlet boundary condition. Recall that the opera- 
tor L with the Dirichlet boundary condition possesses an infinite sequence 
of eigenvalues, all of which have negative real part. By the Krein-Rutman 
theorem, the principal eigenvalue-the eigenvalue with largest real part — 
is real and simple, and the corresponding eigenfunction does not change 
sign The same is true for L, the formal adjoint of L with the Dirichlet 
boundary condition. Furthermore, the principal eigenvalues of L and L co- 
incide. Let (pQ > denote the principal eigenfunction corresponding to 
the principal eigenvalue for L. Normalize it by (pQ {x)dx = 1. Abusing 
notation, we will also let (pQ denote the measure with density (pQ . (The 
measure (p^ is the so-called quasi-invariant distribution for the original dif- 
fusion corresponding to L with killing at the boundary. That is, one has 

IE~5 (/(>'(*)) |rz) > t) = f{x)(p^{x)dx, for all t > 0, where td is the first 

4'o 

exit time of the diffusion Y{t) from D and E~ denotes the expectation for 

the diffusion killed at the boundary and starting from the distribution (pQ.) 
Theorem 1. Consider the operator C in the case that the jump measure is 
given by u = (pQ , where (pQ is the normalized principal eigenfunction for the 
formal adjoint L of L with the Dirichlet boundary condition. Let {A^}J^q 
denote the eigenvalues for L with the Dirichlet boundary condition, labeled 
so that Re Xn+i — > '^'^^ {cpn}^=o be a corresponding sequence 
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of eigenf unctions. Then the eigenvalues for C are and and a 

corresponding sequence of eigenf unctions is given by 1 and {^nlj^i- In 
particular, 

ji{^) = Re (Af). 

Furthermore, (j)^ is the invariant probability measure for the diffusion with 
random jumps from the boundary corresponding to C In fact, (p^ is the 
unique fixed point for the map Inv: V{D) — > V{D) defined in Theorem BP. 

In order to make the spectral analysis tractable when the jump measure v 
is not the special measure considered in Theorem [H we will need to assume 
that the operator L with the Dirichlet boundary condition is self-adjoint, 
although C will still not be self-adjoint, as we now explain. If the first- 
order term b in the operator L is of the form b = aVQ, then the operator 
L can be written in the form L = ^exp(— 2Q)V • aexp(2Q)V. Since we 
can replace Q by Q + c, where c is a constant, without changing L, we will 
assume without loss of generality that J^) exp(2(^)6?x = 1. Let fi^ev 

denote 

the probability measure exp(2Q)dx. In this case, the operator L with the 
Dirichlet boundary condition, considered as an operator on iL^(Z), /ij-cv), is 
symmetric on the domain of smooth functions vanishing at the boundary 
and is self-adjoint on an appropriate domain of definition. (The diffusion 
process in D killed at the boundary, corresponding to L with the Dirichlet 
boundary condition, is reversible and the normalized reversible measure is 
/irev! whence the notation //rev) The operator C, on the other hand, will 
never be self-adjoint. Indeed, a straight forward calculation shows that the 
adjoint operator (with respect to Lebesgue measure) £ of £ is defined on a 
domain which includes {v G C'^{D) n C{D) : w = on dD^, and for such 
functions one has Cv = Lv — {J^ Lv)i>, where L = • aV — bV — V • 6. 
In the case that L is self-adjoint, if one takes the adjoint of C with respect 
to /Urev, then the adjoint is defined on the above class of functions by Cv = 



6 IDDO BEN-ARI AND ROSS G. PINSKY 

We will begin with a key theoretical result, which will be mined to obtain 
more concrete results. Before we can state the theorem, we need some 
additional notation. The eigenvalues of the self adjoint operator L are real 
and negative. We will denote them by {A^};^g, labelled in nonincreasing 
order. Denote the corresponding eigenfunctions by {0^}^O5 normalized by 
jjj (f)!^ djjL^eY = 1) > 0, and (f)^ > 0. Let 

(1.2) Fr,= j cP^dnrev and Gn{l^) = [ (P^du. 

J D J D 

Let {A^I^Q denote the collection of distinct eigenvalues among {\n}'n=o^ 
labelled in decreasing order. We will sometimes need the following assump- 
tion. 

Assumption 1. The Fourier series Yl'^=oTb4>ni.^) converges uniformly 
and absolutely. 

Theorem 2. Assume that the operator L with the Dirichlet boundary con- 
dition is self-adjoint. Let 



n=0 " 

Let dn denote the dimension of the eigenspace corresponding to the n-th 
distinct eigenvalue of L. Assume either that v possesses an L^(D, djirev)- 
density or that Assumption 1 on the operator L holds. 

Then the set of nonzero eigenvalues of C and their multiplicities are given 
as follows: 

i. The set {X : E^{X) = 0} - {A^}~^ ^ consists of simple eigenvalues; 

ii. For each n = 1,2 - •• , the following rule determines whether is an 
eigenvalue, and if so, specifies its multiplicity: 

If dn = 1 and neither F^ = nor Gm{i^) = 0, for the m satisfying = 
A^, then A^ is not an eigenvalue. Otherwise, A^ is an eigenvalue and its 
multiplicity is specified as follows: 

If Gm{i^) 7^ for some m such that = A^ and F^ ^ for some m such 
that A^ = A^, then the multiplicity is dn — I; 
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If Gm{i^) = for all m such that = and Em 7^ for some m such 
that = A^, or if Gm{^) 7^ for some m such that A^ = A^ and Fm = 
for all m such that A^ = , then the multiplicity is dn; 
If Gmii^) = for all m such that A^ = A^ and Fm = for all m such 
that A^ = A^, then the multiplicity is dn if Ey{K^) 7^ and is dn+i if 
E.{K^) = 0. 

Furthermore, even without Assumption 1 or the density condition on u, 
the set of eigenvalues of C includes those listed in (ii). 

Note that the complete characterization of the spectrum in Theorem [2] 
always holds if the jump measure v possesses an L^(Z), /irev)-density. If the 
operator L on D satisfies Assumption 1, then it holds for all jump measures 
y G V{D). The following theorem collects some sufficient conditions for 
Assumption 1 to hold. 

Theorem 3. i. If d = 1, then Assumption 1 holds. 

ii. Let d = 2 and let L = exp(— 2(5)V-a exp(2(5) V satisfy Q = log y^det{a^^) 
(in which case L can be considered as ^Aj\/, where Am is the Laplacian of 
a Riemannian manifold with metric a). Then Assumption 1 holds. 
Hi If d < 3 and the eigenf unctions {(/>^}5^q are uniformly bounded, then 
Assumption 1 holds. 

Remark. A direct calculation (see the proof of Proposition [1] below) shows 
that the eigenfunctions {(1)^}'^=^ are uniformly bounded for L = in 
D = (0, 1)'^; however such a bound does not hold if D is a sphere [2j. 

As a first application of Theorem [2l we identify a class of jump measures 
V for which all the eigenvalues of C are real. The analysis of the spectrum 
in this case turns out to be more tractable. 

Theorem 4. Assume that the conditions of Theorem\^ are in force, and let 
En and Gn{i^) be as in (jl.2p . Assume also that the jump measure v satisifies 
one of the following two conditions: 
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i. FnGn{v) > 0, for all n > 1, or FnGnii') < 0, for all n > 1; 

ii. FnGn 7^ for at most two values of n. 
Then all the eigenvalues of L are real. 

When the nonzero eigenvalue with the largest real part is real, we can 
prove an upper bound on the eigenvalue spectral gap, 71 (z^), of £. 

Theorem 5. // the jump measure v is such that the nonzero eigenvalue of 
C with the largest real part is real, then 

Remark 1. Theorem [S] holds regardless of whether the operator L is self- 
adjoint; however, if L is not self-adjoint then we have no way of determining 
whether the nonzero eigenvalue of C with the largest real part is real. 
Remark 2. As is well known, Ag* gives the exponential rate of decay in t 
of the probability that the diffusion Y{t) in D corresponding to L has not 
yet hit the boundary by time t; that is, lim^^oo 7 log -Px'(''"D > t) = ^0 , 
where td is the first exit time of the diffusion from D. Now since 71(1^) 
gives the exponential rate of convergence of the distribution of the diffusion 
with random jumps to its invariant measure, and since the jump mechanism 
only comes into affect after time td, Theorem [5] might seem (at least at first 
blush) counter-intuitive. 

The normalized reversible measure fj-rev, with respect to which L is self- 
adjoint, plays a distinguished role as the jump measure. In particular, in 
this case the spectral gap can be given by a variational formula. 

Theorem 6. Assume that the operator L with the Dirichlet boundary con- 
dition is self-adjoint. Let {A^},^q denote the eigenvalues of L with the 
Dirichlet boundary condition. Let the jump measure be the normalized re- 
versible measure firev 

1. All the eigenvalues of C are real. 
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11. 



where the infimum is over functions u ^ satisfying u\g£i = Jj^ udfirev 
0. The infimum is attained at a function Umin which satisfies the 
equation Lu = ji{prev)u + C, for some constant C , and the eigen- 
function vi for C corresponding to the eigenvalue jiinrev) is given 



by Vi = Umin + 



c 

71 {^lrev) ' 



111. 



More precisely, consider the function 



E^ir<,v (A) - 



n=0 " 

which is increasing for A G (Af ,A^). // the equation E^^^^{X) = 
possesses a root in (Af,A^), then ji{firev) is equal to this root. 
Otherwise, ji{fJ,rev) = Af . In particular, such a root will exist if 
Fj = (f)^d^rev 7^ 0, for some j G {1, • • • , Uq}, where ko = max{n : 
An = Af}. If Fj = J^(j)fdfj.rev = 0, for all j e {!,••• ,^o}, then 
71 ifJ-rev) > Af if and only if 

Remark. Consider the diffusion process corresponding to L as in Theorem 
O with reflection at the boundary in the conormal direction an, where 
n denotes the inward unit normal to D. The process is reversible and it 
corresponds to a self-adjoint operator on L^(D,;Urcv) which is an extension 
of L with the Neumann boundary condition Vu • an = on dD. For this 
process, /irev is the invariant measure, and the rate of convergence to /irev is 
given by the largest nonzero eigenvalue, A^^. This eigenvalue is given by the 
variational formula in part (ii) of Theorem [6l but with the infimum being 
taken over functions u satisfying ud^^^^ = (without the additional 
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restriction that u\q£) = 0). The infimum is attained at the eigenfunction(s) 
corresponding to Xi , and it is known that any such function does not vanish 
identically on dD. Thus, it follows from part (ii) of Theorem [6] that > 
7i(/t^rev)- Therefore, the rate of convergence to equilibrium is greater for the 
diffusion with random jumps with jump measure //rev than for the reflected 
diffusion, whose invariant measure is /trev 

Here is an application of condition (jl.3p in part (iii) of Theorem [6l 

Proposition 1. Consider the operator in the d- dimensional unit cube, 
D = (0, 1)^, and let the jump measure be Lebesgue measure, l^, on D. One 
has X^ = -^ and Af = 

ilfd< 10, then 71 (Zd) = Af . 
ii. Ifd> 11, then Af < 71 (/d) < A^. 

Remark. Note that 71 (/d) decreases to —00 as d — > 00. Thus, for Brownian 
motion in the d-dimensional cube with random jumps from the boundary 
with normalized Lebesgue measure as the jump measure, the rate of conver- 
gence to equilibrium becomes arbitrarily fast as the dimension increases. 
This is because starting from any point, the distribution of the hitting 
time of the boundary converges to the 5-measure at as d — > 00, which 
means that as d ^ 00, the process constantly gets redistributed according 
to Lebesgue measure after arbitrarily small intervals of time. In contrast to 
this, consider Brownian motion in the d-dimensional unit cube with normal 
reflection at the boundary. The rate of convergence to equilibrium is gov- 
erned by the largest nonzero eigenvalue of the Neumann Laplacian, which 
is Xi = independent of d. Similarly, consider Brownian motion in the 
d-dimensional unit cube, conditioned never to hit the boundary [8]. This 
process corresponds to the /i-transformed operator (^A — X^y^^ . The rate 
of convergence to equilibrium is governed by the largest nonzero eigenvalue, 
which is Af - A^ = = -|7r2, independent of d. 

We have the following result for the one-dimension Laplacian. 
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Proposition 2. Consider the operator in the interval (0,1). 

i. // the jump measure is deterministic; that is, v = 6p, for some p G 
(0, 1), then 71(1/) = Af = -2Tr^; 

ii. // the jump measure v is such that the nonzero eigenvalue of C with 
the largest real part is real, then 

7i(z.) =Af = -27r2. 

Remark. Part (i) above was shown in [3] and [5]. (Actually, — ^ was 
obtained in [3], because a certain cancellation was not taken into account. 
The correct result appears in [4J.) Our proof is completely different. A direct 
calculation shows that Fi = 0; thus, by Theorem [21 71 (z^) > Af = — 27r^, 
for all z^. In a preprint version of this paper, we made the conjecture that 
71 (i^) = Af = — 27r^ for all u. This conjecture has now been established by 
combining part (ii) with a very recent result [7] which states that in the case 
of on an interval, all of the eigenvalues are real, for all jump measures 
u. 

The next result shows that for the Laplacian on a square in R^, one can 
find a deterministic jump measure for which 71(1^) 7^ Af. 

Proposition 3. Consider the operator in the square (0, 1)^. Then there 
exists a jump measure v of the form bx^^ for some € (0, 1)^, for which 
7(z^)>Af = -^. 

Remark. Consider L = in the cube (0, 1)*^, d > 1. The proof of Propo- 
sition [T] shows that -Fi = 0; thus, by Theorem [21 7i(i^) > Af . Combining 
Proposition [3] and Proposition [H it follows that for d > 11 or d = 2, there 
exists a jump measure u for which 71(1^) > Af . Presumably, this holds for 
all d > 2. Conversely, by Theorem [1] it follows that for all d > 1 it is also 
always possible to find a u for which 71 (z^) = Af . 

In all of the examples given so far, 71 (z^) > Af . The following result 
shows that such is not always the case. 
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Proposition 4. Let ko = max{n : = Af}. If Fj ^ 0, for some j G 
{1, 2, • • • , A;o}, then it is always possible to find a jump measure v for which 
7i(i/) > Af and it is always possible to find a jump measure v for which 
< Af . 

Remark. One can check that L = on (0, 1), where 5 is a nonzero 

constant, is an example where Proposition U] is apphcable. 

We conclude this section with several open questions. 
Question 1. In a preprint version of this paper, we asked whether all the 
eigenvalues of C are real in the case that L is self-adjoint. It the very 
recent paper [7] it was shown that for L = in a ball in , there exist 
deterministic jump measures, that is measures of the form v = 5x, for which 
some of the eigenvalues are not real. However, these non-real eigenvalues 
do not have maximal real part. We still ask whether the eigenvalue with 
largest real part is real in the case that L is self-adjoint. 
Question 2. Does the inequality 71 (z^) < Xq hold for all v when L is self- 
adjoint? What about for general L? 

Question 3. Does a lower bound exist for 71(1^) in terms of the eigenvalues 
{A^}-o ofL? 

Question 4. What can be said about the continuity properties of 71(1^) as v 
varies over V{D)^ the space of probability measures on D with the topology 
of weak convergence? 

Remark. Note that if Question 1 is answered affirmatively, then Theorem 
[5] shows that the answer to Question 2 is affirmative in the case that L is 
self-adjoint. 

The proofs of the results stated in this section are grouped thematically 
and proved in the sections that follow. 

2. Proof of Theorem [T] 

We first prove the statement concerning the eigenvalues and eigenfunc- 
tions. Let 0^, n > 1, denote an eigenfunction for L with the Dirichlet 
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boundary condition, corresponding to the eigenvalue A^. Integration by 
parts yields 

Jd Jd Jd Jd 

from which it follows that J^(f)^4>Qdx = 0. Since 4>n\dD = 0, it follows 
that (p^ is in the domain of C, and we conclude that is an eigenvalue 
for C. Of course the function 1 is an eigenfunction for C corresponding to 
the eigenvalue 0. Thus, to conclude the proof, it is enough to show that 
the equation Lip = Xip with ipcp^dx = V'bz) = c 7^ 0, has a solution only 
when A = 0. Let V be a solution to the above equation. Let da denote 
the Lebesgue surface measure on dD. Integrating by parts and using the 
boundary condition, we have 

\D f Id,/, _ / j,taD 



X^c = X^ / (t)^il)dx= / ipLcp^dx 
= I (pQ Lipdx + I -ipaV<pQ ■ nda — f ipfp^b ■ nda 

Jd JdD 2 JdD 

= X (p^ijjdx + c / L(j)^dx = cX + cXq, 
Jd Jd 

where n is the unit outward normal of D at dD. Therefore A = 0. 

We now turn to the statement concerning the invariant measure. We 
denote by the operator from to given by {G^ f){x) = G^{x,f). 
We denote its adjoint on by G^. We have {G^g){y) = G^{g,y). We 
now prove that G^ is compact. For e > 0, let denote the operator on 
defined by 

{KJ){y) = p''isj,y)ds. 

An argument similar but simpler than the one given in the proof of Lemma 
1 in [1], based on the continuity of on [e, e^^] x D x D, shows that is 
compact. Now, 

\\GDf-KJ\\,= / I / p''is,f,y)ds+ / p''{.sj,y)ds\dy 
Jd Jo Jt~^ 

(2.1) <e||/||i + ||/||i / supPf (n >s)(is. 

Je-i x&D 
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By O Theorem 3.6.1], limg^oo ^^ogsup^^D Px{ti > s) = \q . Thus, there 
exists a constant C > such that for ah e sufficiently small. 



-1 x£D 

Therefore, it follows from (12. ID that 



supPf (ri > s)ds < Ce^< 



^ 

\\GD - K,\\i <e + Ce^'~' ^^0. 

Consequently, is compact. 

Assume now that m G P is a fixed point for Inv. Since dm{y) = 
G (m,y) , _ G m fg^jQ-^yg ^j-^at m has density in . Therefore we 

G-D(m,l) ^ G°(m,l)' 

may consider m as an eigenfunction for G-^, corresponding to the eigenvalue 
A = {m,\). Let (j)^ > denote the principal eigenfunction for L with 
the Dirichlet boundary condition corresponding to the eigenvalue A^. Since 
G^ = (-L)~i, we have G^0o = -(A^)-^?- Therefore, 

A / m(t)^dx = [ G^m(t)^dx = [ mG°<i)^dx = -{XqY^ f mtp^dx. 
Jd Jd Jd Jd 

Since m and (pQ are nonnegative, it follows that A = — (A^)^^. Since 
G^ is compact, the Krein-Rutman theorem guarantees that the eigenvalue 
— (A^)~^ for G^ is simple. Since (po and m are both eigenfunctions for G^ 
corresponding to the eigenvalue — (Ao')~^, and since mdx = (podx = 1, 
we conclude that m = (bn. 



□ 



3. Proofs of Theorem [2], Theorem m Theorem [6], Proposition □ 

AND Proposition HI 



Proof of Theorem [2l A number A € C — {0} will be an eigenvalue if and 
only there exists a function v satisfying Lv = Xv and v\gD = Jj^vdv. Let 
u = V — c, where c = v\qd- Then u satisfies Lu = Xu + K, where K = Ac. 
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On L'^{D 

)/^rev)) the function u can be represented in the form 

oo 

(3.1) u=Y,CnC^n, 

n=0 

for unknown constants {Cn}'^=Q-, and the constant function 1 can be repre- 
sented by 

oo 

(3.2) 1 = ^F„0^, 

n=0 

where {-F„.}^o ^ ™- (112]) • Since u is a smooth function vanishing on 
dD, it is in the domain of the self-adjoint operator L acting on L^{D, firev)- 
Thus, from (13.11). it follows that 



(3.3) Ln = ^A^C„0f. 

n=0 

From (j3.ip - ()3.3p along with the fact that Lu = Xu + K, we obtain 

(3.4) CnX^ = \Cn + KFn, n>G. 

We first show that the condition E^{X) = is necessary and sufficient 
for A {^n In^o eigenvalue. Since we are now assuming that A is 

not in the spectrum of L, we may assume that K ^ 0. Indeed, if K were 
equal to 0, then v would vanish on dD and consequently it would be an 
eigenfunction for L. This would mean that A = A^, for some n. From p.4p 
we obtain 



rf — " 



AP-A' 



and conclude that 
(3.5) 



n=0 " 



In order that v be an eigenfunction, v must satisfy v\qd = J^v du = c. 
Since u = v — c, we require that Jj^ u dv = 0. If has an L^i^D, ^rev 

)-density, 



16 
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then using p.Sp and taking inner products shows that jj^u dv = if and 
only if 



where G'n(i^) is as in (11. 2p . Alternatively, if Assumption 1 holds, then the 
formula for u in (13. 5p holds not only in L'^{D , firev), but also pointwise, and 
from the bounded convergence theorem it follows again that j^u dv = Q \i 
and only if (13. 6p holds. We have thus shown that the condition Ey{\) = 
is necessary and sufficient for a nonzero A {A^},^]^ to be an eigenvalue. 
Furthermore, as the method uniquely specifies the corresponding eigenfunc- 
tion (up to a multiplicative constant), it follows that the multiplicity of such 
an eigenvalue is 1. 

We now consider the possibility that A = A^^^ is an eigenvalue, where 
no is a nonnegative integer. Let 5„o denote the -dimensional eigenspace 
corresponding to the eigenvalue A^^^ of L. Let S'^^{v) = {u; G Sn^ '■ J^) wdv = 
0} and let S^^ = G 5„ : f^wdfi^cv = 0}. Clearly, each of these latter 
two spaces is either (d^g — l)-dimensional or dng-dimensional. 

Consider first the case that S^^ is {dno — l)-dimensional. There exists 
an mo such that A^^^ = A^^ and = (pmod^J-rev 0. But then (j3.4p 
will hold with n = mo and A = An^ if and only if X = 0. But if X = 0, 
then u = V, v\g£) = vdu = and Av = A^^f . Thus, v belongs to S^^{fi). 
Consequently, the multiplicity of A^^ will be either dn^ — 1 or dn^ , depending 
on which of these numbers is the dimension of S^^{fj,). In particular, if 
no = 0, then dn^ = 1 and S^^ = S^^{u) = {0} since (pQ > 0. Thus, 
Aq* = Aq* can never be an eigenvalue. 

Now consider the case that S^^ is d„g-dimensional. In this case, Fm = 0, 
for all m such that A^ = A?^. We first look for eigenfunctions for which 



(3.6) 
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K ^0. Solving (fO]) gives 
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Cn = ^3a' ^ ^^^^ ^^^^ 7^ ^no! 

no- 



C„ is arbitrary, for all n such that = 



Writing C„ = Kcn, for n such that = A^^, and employing the same 
reasoning as in (|3.5p and (|3.6p yields 



(3-7) E E C"G"(-) = 0- 







There are two cases to consider — when S^^{i') is (d„,) — l)-dimensional and 
when it is d„Q -dimensional. In the latter case, Gn{i^) = 0, for all n satisfying 
= A^^. Thus, (j3.7p reduces to E^{A^^) = 0. If this equation is satisfied, 
we obtain one eigenfunction with K ^ 0, and if it is not satisfied, we obtain 
no such eigenfunctions. Since S^^{i') is -dimensional, there are also dn^ 
additional linearly independent eigenfunctions with K = 0. Thus, the mul- 
tiplicity is either d„g -|- 1 or d„,), depending on whether or not E^{A^^) = 0. 

Now consider the case that S^^{v) is — l)-dimensional. Since we 
may choose the orthonormal basis {0m}{m:AO=AO } corresponding to the 
eigenspace 5„q however we like, we may assume without loss of generality, 
that Gm{^) = 4>mdv = 0, for all but one of the m for which = A^^. 
Denote the single m for which this is not true by mo- Then (j3.7p reduces to 



EFnGn{v) „ ^ ^ _ n 



no 



The above equation is uniquely solvable for , and thus yields one eigen- 
function with K 0. Since S^^{u) is (d„o — l)-dimensional, there are also 
— 1 additional linearly independent eigenfunctions with K = 0; thus the 
multiplicity is dnQ- □ 
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Proof of Theorem [4l By Theorem [21 a complex number A = a + with 
P will be an eigenvalue for C if and only if 

, A^ - A " 

n=0 " 



We can rewrite this as 



(3.8) 



n=0 

FnGn{v)f3 



(A^-a)2 + /32 
^ FnGn{v)(3 ^ 



Clearly, the two equations in (|3.8p hold if and only if the following two 
equations hold: 



(3.9) 



E FnGn{v) _ „ 

(A^ - af + /32 - 



Since FoGo(j^) is always positive, neither equation in (j3.9p can hold if FnGn{i^) > 
0, for all n > 1. Consider now either the case that FnGnii^) < 0, for all 
71 > 1, or alternatively, the case that FnGn{i^) is nonzero for no more than 
two values of n. Rewriting (j3.9p as 

F^Goiu) ^ FnGnjv) Ag ^ 
(Ao^ - a)2 + /?2 + (A^ - a)2 + /32 ' 

FoGo(r/) F„G„(t/) _ 

(Ao^-a)2 + /?2+A.(A^-a)2 + /32 ' 

it follows that the two equations in (|3.9p cannot hold simultaneously. □ 
Proof of Theorem [6l i. Since Gni^J^ev) = Fn, it follows from part (i) of 
Theorem m that all the eigenvalues of C are real. 

ii. By part (i), 71 = 7i(/Urev) is itself an eigenvalue; let (pi denote a corre- 
sponding eigenfunction. Let ^pl = (pi — c, where c = (pi\dD = Jd 4'idfJ-Tev 
Then ^pl\QD = /^'^id/^rev = and Lipi = ^iipi + 71c. Multiplying this 



ELLIPTIC OPERATORS WITH NONLOCAL BOUNDARY CONDITION 19 

equation by tpi exp{2Q) and integrating by parts gives 

rev 

Jd V'l^^^rev 

On the other hand, consider the quotient 2 /D(V«aVM)dMrev ^ standard 
methods, the infimum of this quotient over functions 7^ n € Hq{D) sat- 
isfying u\d£) = f^udfj^rev = exists. We denote this infimum by —T > 0. 
To identify the minimum, we use a Lagrange multipher and vary the quan- 
tity ^ Jjj{\7uaVu)dfircv + k Jj^u^dfirev over functions u satisfying the above 
restriction, where A; is a free parameter. A minimizer ip must satisfy the 
equation q{Lip — kip) d^rcv = 0, for all q satisfying the above restriction. 
From this one concludes that Lip = kip + C, for some constant C. Multi- 
plying this equation by ip, integrating both sides with respect to dfi^ev, and 
integrating by parts, one finds that k = T. Letting (p = ip + jr, '^t follows 
that (p satisfies L(p = Tcp and (p\dD = J^, </'d^rev 

in. By part (i) and the definition of 7i(^rev), it follows that 7i(/irev) is the 
largest nonzero eigenvalue of C. And then by Theorem [5] it follows that 
liifJ-rev) < A^. In Theorem [21 note that when then GnifJ-rev) = 

Fn. Consequently -E„^^„(A) = Xl^o 1^31- Since -E'^^^^(A) is continuous for 
A G (Af , A^), since ^^„„((A^)-) = 00 and since ^^^.^((Af )+) = -00 holds 
if Fj = fjj<pfdfirev / 0, for some j E {1,2, ■•• ,ko}, it follows that E^j^^^^ 
possesses a root in (Af, A^) if Fj / for some j G {1, 2, • • • , ko}. It now 
follows from Theorem [2] that Af < '~^i{iirev) < A^, with strict inequality if 
Fj 7^ 0, for some j G {1, 2, • • • , k^}. Furthermore, since E^^^^ is increasing 
on (Af , A^), if follows that in the case that Fj = for all j G {1, 2, • • • , k^}, 
the strict inequality will hold if and only if E^j_^^^{\^) < 0. This inequality 
can be rewritten as p.3|) . □ 

Proof of Proposition [H By Theorem [61-iii, 7i(/iret;) < Ag*. To prove the 
rest of the proposition, we apply (|1.3p from Theorem [6l The complete, or- 
thonormal sequence of eigenfunctions on L'^(D,l^) for iA on D = (0,1)'^ 
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with the Dirichlet boundary condition is given by 
{22 J^^^j^ sinnjvrxjlj^ ... The corresponding eigenvalues are 
{— ^ X^^^]^ n|}5^ ... j^^^;^. We will denote these eigenfunctions and eigen- 
values respectively by ... „ , and A£ ... „ ,. We have 



•ni,--- ,nd' 



(3.10) Fn„ 



D 



dx 



23" 



TT^n 



if rij is odd for all j; 



0, otherwise. 



In the present context, the terms Fq, and Af appearing in Theorem[6]are 
given respectively by = Ai,...,i = and A„i,...,„^ = _l£^±|)z:_^ 

where (ni • • • , n^) satisfies Yl'i=i — d+l. From (jS.lOp . we have -F„j^... = 



Id ^ni,- ,nd d'fJ-vev = 0, if Yl'j=i nj = d + 1. Thus, (|1.3p is apphcable. 

Using {Arti,...,nd} and {F„^,...,„^} in place of the labeling {A^} and {F„} 
in the inequality (|1.3|) . we find that after cancellations the inequality can be 
written as 



(3.11) 



1 



ni,- ,nd odd 11^=1 n^j (Ej=i n]-d-3 



1 

>3- 



Thus, by (|1.3p . (|3.1ip is a necessary and sufficient condition in order that 
7i(/^rev) > Af, and if the condition does not hold, then 7i(/Urov) = Af. 
Denote the left hand side of p. lip by Hd- If one considers H^+i, but restricts 
the summation to those multi-indices (ni,-- - for which n^+i = 1, 

the resulting quantity is Hd. Thus the left hand side of (j3.1ip is monotone 
increasing in d. A direct calculation shows that the inequality in (j3.1ip does 
not hold if d = 1. On the other hand, by considering the contribution to the 
left hand side of (j3.1ip only from those multi-indices satisfying Ej=i nj = 
d + 2, it is easy to check that the inequality in (j3.1ip holds if d > 15. 
From these observations we conclude that there exists ad* G [2, 15] such 
that 7i(/irev) > Af, if d > d*, and 7i(/Urov) = Af, if d < d*. Numerical 
calculations shows that in fact d* = 11. □ 
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Proof of Proposition |4l Without loss of generality, assume that Fi > 
0. Choose with density iy^{x) = c±{4'q{x) it e(/)f (x)), where e > is 
sufficiently small so that z^^(x) > 0, for all x £ D, and where c-t > is a 
normalizing constant so that z/^ is a probability density. (It is possible to 
choose such an e > because the Hopf maximum principle guarantees that 
the normal derivative V4>q{x) ■ n ^ 0, for x G dD.) Recall that G„(z^^) = 
(/>^ di^^. Thus, Go{u^) = c±, Gi{i^^) = ±c±e and G„(z^=^) = 0, for n>2. 
From the definition of Ej^± in Theorem [2l we have E^±(X) = Yn~T ^ 'ur~k- 

Aq - — A Aj^ — A 

Thus, E^-{\) for A e (Af ,A^). By Theorems d and [5l 71(1/") < A^. 
Thus, we conclude from Theorem [2] that 71(1/^) < Af. Since E^+{{Xq)~) = 
00 and £'^+((Af)^) = —00, £'i,+ (A) possesses a root in (Af, A^). Thus, by 
Theorem [21 71(1/+) > Af . □ 



By assumption, 71(1^) is a real eigenvalue for C. We need to show 
that 71(1^) < A^. Let u denote a corresponding eigenfunction, and let 
c = u\d = udv. We first show that 71 (z^) 7^ A^. Assume to the contrary. 
In this case, c 7^ 0. Indeed, otherwise (pQ and u would both be eigenfunctions 
for the principal eigenvalue A^ of L with the Dirichlet boundary condition. 
Furthermore, and u would be linearly independent since (p^ does not 
change sign, whereas udv = c = 0. This would then contradict the sim- 
plicity of the principal eigenvalue A^. Integrating by parts twice, exploiting 
the form of the reversible operator L and the reversible measure, we have 



4. Proof of Theorem [5] 




'rev — 





'rev 1 



which reduces to 



(4.1) 
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However, 



(4.2) 



/ uaV(j)Q ■ nexp{2Q)da = c aV(j)Q ■ nexp{2Q)da 
JdD JdD 



IdD JdD 

rev ■ 



' D J D 

Now (j4.ip and (j4.2p give c = 0, which is a contradiction. 

Now we show that 71(1^) ^ A^. Assume to the contrary. By the Feynman- 
Kac formula, u{Y{t A r^))) exp(— 7i(i/)(t A td)) is a martingale. Thus, 

(4.3) Efn(y(t Ar^))exp(-7i(z.)(t ATfl)) =n(x). 

Since 7i(j^) > A^, we have exp(— 7i(z^)r£)) < 00 [9, chapter 3]. Thus 
letting t ^ 00 in ()4.3p and applying the dominated convergence theorem 
gives 

(4.4) ti(x) = cEf exp(-7i(i/)TD). 

It follows from (j4.4p that c 7^ 0. Integrating both sides of (|4.4p against 
now gives 

(4.5) E^exp(-7i(z.)ri5) = l, 

which is a contradiction. □ 

Remark. If one does not assume that the nonzero eigenvalue with largest 
real part is real, the calculation in the above proof can be made with ^i{v) 
replaced by Ai(i^), where Ai(i^) is an eigenvalue for C whose real part is 
7i(i^). One arrives at (j4.5p with 71(1^) replaced by Ai(i^). However, since 
Ai(z^) can be complex- valued, (j4.5p no longer constitutes a contradiction. 

5. Proof of Proposition [2] 
i. The eigenvalue problem for C is 



(5.1) 



\u" = Xu in (0,1); 
u{p) = u{0) = n(l). 
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Every solution u to the differential equation in ()5.ip is given by 
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u{x) = A cos Kx + B sin nx, 

where A = — ^k^. In order that such a solution also satisfy the boundary 
condition in (jS.ip . the following system of linear equations must have a 
nontrivial solution: 

j4 (1 — cos Kp) —B sin Kp = 0; 
(5.2) A(l - cosk) -5sinK = 0. 

The determinant of the linear system above is 
(5.3) 

— (1 — cos Kp) sin K + sin Kp{l — cos k) 

= sin k{1 — p) — sin k + sin np 

„ . k(1 — n) k(1 — p) k(p+1) . 

= 2 sin ^ cos — + 2 cos — sin k p - 12 

2 2 2 ^ 

k(1 -p) / k(1 -p) k(p+ 1) 
= 2 sm cos cos 



^ . k(1 — p) . K . Kp 

= 4 sm sm — sm — . 

2 2 2 

Since sinx = if and only if x = nn for some integer n, the solutions k of 
(j5.3p are all real and are given by 

27rn 27rn 

, 27rn, , n G Z. 



17 7 
— P P 

Therefore, the eigenvalues for C are 

-2 , -27rV, n G N. 

(1 — pj^ p^ 

Thus, the non-zero eigenvalue with maximal real part is —2tt'^. 

a. By assumption, the nonzero eigenvalue with largest real part is real, and 

we know that it is negative. Thus, 71 (i^) is the largest negative number 7 
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for which there is a solution to the fohowing problem: 



(5.4) 



( 

lu" = -fu in (0,1); 



u{0) = u{l) = Jq u du. 



Every solution to the differential equation in (j5.4p with 7 < is of the form 

u{x) = A cos Kx + B sin nx, 

where 7 = for some k ^ R — {0}. In order that (j5.4p have a solution, 

the following system of linear equations must have a nontrivial solution: 

— J cos KX du^ —B j sinKx dv = 0; 
(5.5) A(l - cosk) -SsinK = 0. 

The determinant of the linear system above is 

— (1 — / cos dz^) sin K + (1 — cos k) / sninx dv 
Jo Jo 
1 .1 
sin«;(l — x) dz^ — siuK + / sinux dv 
Jo 

. k{1-x) k{1-x) k{x + 1) . k{x - 1) 
2 / sm cos dv + z 1 cos sm dv 



^ ^ JO 

. k(1-x)/ k(1-x) k(x + 1)\ 
2 / sm cos cos dv 



k(1 — x) k kx 

= 4 / sin sin — sin — dv. 

Jo 2 2 2 

Note that sin ^^^i^"^ sin f sin ^ > 0, for k G (0, 27r) and x G (0,1), while 
the reverse inequality holds for k € (— 27r,0) and x £ (0,1). Thus, it fol- 
lows that (j5.5p has no solution k € (— 27r,27r) — {0}. On the other hand. 



sm 



sin f sin = 0, for k = ±2tt. Thus, it follows that 7 = -\{2TTf 



27r is the largest real nonzero solution to ()5.4p . □. 
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6. Proof of Proposition [3] 

By Theorem [3] (part (ii) or part(iii)), Theorem [2] holds for all jump 
measures i'. We will show that for an appropriate v one has Ey[\) = 0, for 
some A > Af. 

We use the notation and the calculations in the proof of Proposition [TJ 

Let = (^(i,!)- We have G„,,„2(i/o) = Jo^t^n^^n^ d^o = 'Pni,n2ihh) = 
2 sin(^7r) sin(^7r). Then from the definition of Ey^{X) it follows that 



C 



(6-1) ^ ^ sin(^^7r)sin(^^7r) 



mi ,m2- 



(2mi + l)(2m2 + 1) ((2mi + 1)2 + (2^2 + 1)2 + ^\) 



for an appropriate negative constant C . We will show that the equation 
Ey^,{\) = has a root A G (Af ,A^) = (-fvr^ -vr^). Note that ^^^((-Tr^)-) = 
oo. Thus, it suffices to show that ^) < 0. This can be checked using 

a program such as Mathematica, or alternatively, by a page and a half of 
estimates which we refrain from reproducing here. 



7. Proof of Theorem [3] 

For all three parts of the theorem, we will need the following comparison 
result. By the mini- max principle [lOj, one can compare the eigenvalues 
{A^IJ^Q of L in D to those of in (0, l)'^, and conclude that there exist 
ci,C2 > (depending on L and D) such that ciA„ < A^ < C2A„, where 
{AnjJ^o ^'^^ eigenvalues for ^ A in (0, l)'^, labelled in nonincreasing order. 
i. It is known that the eigenfunctions {(/>^}J^q are uniformly bounded [6l 
pp. 270-273]. Thus, using the Cauchy-Schwarz inequality, it is enough to 
show that X^n^o (A^p" ^ comparison principle above, it suffices 

to show the above inequality in the case that L = \n D = (0, 1). In 
this case, A^ = _(l!±^.^_ 
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ii. When L is the Laplacian on a Riemannian manifold, it is known that 
|<An I ^ C'jA^l^r, for some C > [2j. Using this and applying the Cauchy- 
Schwarz inequality, it follows that Y1'^=q T^^'ni^)^ converges uniformly and 
absolutely when d = 2 if Yl'n=o l^n converges. By the mini-max principle 
above, it suffices to show the above inequality in the case that L = 
on (0,1)^. But this then follows from Weyl's asymptotic distribution of 
eigenvalues [IQ] which gives Xyi ~ CTi. 

in. As in part (i), it is enough to show that X]jS=o ^ ^^'^ 

comparison principle above, it suffices to show the above inequality in the 
case that L = on (0, 1)^, d < 3. But this then follows from Weyl's 

2 

asymptotic distribution of eigenvalues [TO] , which gives Xn ~ cn ^ , 

□ 
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